CURVATURE FLOWS IN SEMI-RIEMANNIAN MANIFOLDS 



CLAUS GERHARDT 



Abstract. We prove that the limit hypersurfaces of converging curva- 
ture flows are stable, if the initial velocity has a weak sign, and give a 
survey of the existence and regularity results. 
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1. Introduction 

In this paper we want to give a survey of the existence and regularity results 
for extrinsic curvature flows in semi-Ricmannian manifolds, i.e., Riemannian or 
Lorentzian ambient spaces, with an emphasis on flows in Lorentzian spaces. In 
order to treat both cases simultaneously terminology like spacelike, timelike, 
etc., that only makes sense in a Lorentzian setting should be ignored in the 
Riemannian case. 

The general stability result for the limit hypersurfaces of converging curva- 
ture flows in Section [5] is new. The regularity result in Theorem l6.5l — especially 
the time independent C m+2 ' Q -estimates — for converging curvature flows that 
are graphs is interesting too. 
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2. Notations and preliminary results 

The main objective of this section is to state the equations of Gaufi, Co- 
dazzi, and Weingarten for hypersurfaces. In view of the subtle but important 
difference that is to be seen in the Gaufi equation depending on the nature 
of the ambient space — Riemannian or Lorcntzian — , which we already men- 
tioned in the introduction, we shall formulate the governing equations of a 
hypersurface M in a semi-Riemannian (n+l)-dimensional space N, which is 
either Riemannian or Lorcntzian. Geometric quantities in N will be denoted 
by (fla/3)) (Rap-ys), etc., and those in M by (gij), (Rijkl), etc. Greek indices 
range from to n and Latin from 1 to n; the summation convention is always 
used. Generic coordinate systems in N resp. M will be denoted by (x a ) resp. 
(£'). Covariant differentiation will simply be indicated by indices, only in case 
of possible ambiguity they will be preceded by a semicolon, i.e. for a function 
u in N, (u a ) will be the gradient and (u a p) the Hessian, but e.g., the covariant 
derivative of the curvature tensor will be abbreviated by R a /3jS;e- We also point 
out that 

(2.1) Raf3y6;i = Ra0j6;eXi 

with obvious generalizations to other quantities. 

Let M be a spacelike hypersurface, i.e. the induced metric is Riemannian, 
with a differentiable normal v. We define the signature of v, a = cr(v), by 

(2.2) a = g aP v a v^ = (u,u). 

In case N is Lorentzian, a = — 1, and v is timclike. 

In local coordinates, (x a ) and (£*), the geometric quantities of the spacelike 
hypersurface M are connected through the following equations 

(2.3) x% = -o-h t3 v a 

the so-called Gaufi formula. Here, and also in the sequel, a covariant derivative 
is always a full tensor, i.e. 

(2.4) — Xjj — r^x k + r^x^x'J. 

The comma indicates ordinary partial derivatives. 

In this implicit definition the second fundamental form (hij ) is taken with 
respect to — ov. 

The second equation is the Weingarten equation 

(2.5) v? = h\xl, 

where we remember that i/f is a full tensor. 
Finally, we have the Codazzi equation 

(2.6) h ij;k ~ hik-j = RaW a x?x]x 5 k 
and the Gaufi equation 

(2.7) Rijki = crihikhji - huhjk} + Rap^sxfx^ x 7 k xf . 
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Here, the signature of v comes into play. 

2.1. Definition, (i) Let F e C°(f) n C 2 ' a (r) be a strictly monotone cur- 
vature function, where r C R™ is a convex, open, symmetric cone containing 
the positive cone, such that 

(2.8) F \or=° A F ]r >0. 

Let N be semi-Riemannian. A spacelike, orientablcQ hypersurface M C N 
is called admissible, if its principal curvatures with respect to a chosen normal 
lie in r. This definition also applies to subsets of M. 

(ii) Let M be an admissible hypersurface and / a function defined in a 
neighbourhood of M. M is said to be an upper barrier for the pair (F, /), if 

(2-9) F ]M >f 

(iii) Similarly, a spacelike, orientable hypersurface M is called a lower barrier 
for the pair (F, /), if at the points S c M, where M is admissible, there holds 

(2.10) F ]E < f. 
S may be empty. 

(iv) If we consider the mean curvature function, F = H, then we suppose F 
to be defined in R" and any spacelike, orientable hypersurface is admissible. 

One of the assumptions that are used when proving a priori estimates is that 
there exists a strictly convex function x € C 2 (Cl) in a given domain Q. We 
shall state sufficient geometric conditions guaranteeing the existence of such a 
function. The lemma below will be valid in Lorentzian as well as Riemannian 
manifolds, but we formulate and prove it only for the Lorentzian case. 

2.2. Lemma. Let N be globally hyperbolic, Sq a Cauchy hypersurface, (x a ) 
a special coordinate system associated with Sq, and C N be compact. Then, 
there exists a strictly convex function \ £ C 2 (Cl) provided the level hypersur- 
faces {x° — const} that intersect Cl are strictly convex. 

Proof. For greater clarity set t — x°, i.e., t is a globally defined time function. 
Let x — x(£) be a local representation for {t = const}, and U,tij be the 
covariant derivatives of t with respect to the induced metric, and t a , t a p be the 
covariant derivatives in N, then 

(2.11) = Uj =t af 3xTx? +t a xfj, 
and therefore, 

(2.12) t a j3X^ Xj — t a x^j — hijt a v . 

Here, (y a ) is past directed, i.e., the right-hand side in (|2.12p is positive definite 
in fi, since (t a ) is also past directed. 



A hypersurface is said to be orientable, if it has a continuous normal field. 



4 



CLAUS GERHARDT 



Choose A > and define \ — eM > so that 

(2.13) Xa p = AV'M/3 + \e xt t a p. 

Let p G H, be arbitrary, S = {t = t(p)} be the level hypersurface through p, 
and (r] a ) G T P (N). Then, we conclude 

(2.14) e- xt x a p V a rf = \ 2 \if\ 2 + A/,,//'//*' + 2At 0j W, 

where now represents the left-hand side in (|2.12|) . and we infer further 
e^Xa/^V > |A 2 |r?| o2 + [Ae - c^rfrf 



> IX{-\ V °\ 2 + 



(2.15) 

for some e > 0, and where A is supposed to be large. Therefore, we have in 
(2.16) X Q /3>cg Q;3 , c> 0, 

i.e., x is strictly convex. □ 

3. Evolution equations for some geometric quantities 

Curvature flows are used for different purposes, they can be merely vehicles 
to approximate a stationary solution, in which case the flow is driven not 
only by a curvature function but also by the corresponding right-hand side, 
an external force, if you like, or the flow is a pure curvature flow driven only 
by a curvature function, and it is used to analyze the topology of the initial 
hypersurface, if the ambient space is Riemannian, or the singularities of the 
ambient space, in the Lorentzian case. 

In this section we are treating very general curvature flows0 in a semi- 
Riemannian manifold N = N n+ , though we only have the Riemannian or 
Lorentzian case in mind, such that the flow can be either a pure curvature flow 
or may also be driven by an external force. The nature of the ambient space, 
i.e., the signature of its metric, is expressed by a parameter a = ±1, such that 
a = 1 corresponds to the Riemannian and a = — 1 the Lorentzian case. The 
parameter a can also be viewed as the signature of the normal of the spacelike 
hypersurfaces, namely, 

(3.1) <r={v,v). 

Properties like spacelike, achronal, etc., however, only make sense, when TV 
is Lorentzian and should be ignored otherwise. 

We consider a strictly monotone, symmetric, and concave curvature F G 
C 4,a (r), homogeneous of degree 1, a function < / G C 4, "(]?), where fl C iV 
is an open set, and a real function # G C 4 '"(IR + ) satisfying 

(3.2) <P > and 3 < 0. 
For notational reasons, let us abbreviate 

(3.3) /=*(/)• 



2 We emphasize that we are only considering flows driven by the extrinsic curvature not 
by the intrinsic curvature. 
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Important examples of functions 4> are 

(3.4) <P(r) = r, <P(r) = logr, <P{r) = -r -1 
or 

(3.5) $(r)=r*, $(r) = -r~i, k > 1. 

3.1. Remark. The latter choices are necessary, if the curvature function F 
is not homogeneous of degree 1 but of degree k, like the symmetric polynomials 

1 /k 

Hk- In this case we would sometimes like to define F = Hk and not H k , since 

dF 

(3.6) F 11 = 



dhij 



is then divergence free, if the ambient space is a spaceform, cf. Lemma 15.81 on 
page [Ml though on the other hand we need a concave operator for technical 
reasons, hence we have to take the k-th root. 

The curvature flow is given by the evolution problem 

x = —a{$ - f)v, 

(3.7) V J) ' 

x(0) = x , 

where Xq is an embedding of an initial compact, spacelike hypersurface Mo C fl 
of class C 6,Q , <P = 'P(F), and F is evaluated at the principal curvatures of the 
flow hypersurfaces M(t), or, equivalently, we may assume that F depends on 
the second fundamental form (hij) and the metric (gij) of M(t); x{i) is the 
embedding of M(t) and a the signature of the normal v = v(t), which is 
identical to the normal used in the Gaussian formula (|2.3p on paged] 

The initial hypersurface should be admissible, i.e., its principal curvatures 
should belong to the convex, symmetric cone r C W 1 . 

This is a parabolic problem, so short-time existence is guaranteed, cf. [181 
Chapter 2.5] 

There will be a slight ambiguity in the terminology, since we shall call the 
evolution parameter time, but this lapse shouldn't cause any misunderstand- 
ings, if the ambient space is Lorentzian. 

At the moment we consider a sufficiently smooth solution of the initial value 
problem (|3.7p and want to show how the metric, the second fundamental form, 
and the normal vector of the hypersurfaces M(t) evolve. All time derivatives 
are total derivatives, i.e., covariant derivatives of tensor fields defined over the 
curve x(t), cf. 17, Chapter 11.5]; t is the flow parameter, also referred to 
as time, and (£ l ) are local coordinates of the initial embedding xq — xq(£) 
which will also serve as coordinates for the the flow hypersurfaces M(t). The 
coordinates in N will be labelled {x a ), < a < n. 

3.2. Lemma (Evolution of the metric). The metric of M(t) satisfies the 
evolution equation 

(3.8) 9ij = -M*-f)hij- 
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Proof. Differentiating 

(3.9) Qij = (xi : Xj) 
covariantly with respect to t yields 

Qij — \Xi,Xj) \Xi,Xj) 

(3.10) : 

= -2a($ - f){x t , Vj ) = -2tr(# - f)h tj , 
in view of the Codazzi equations. □ 

3.3. Lemma (Evolution of the normal). The normal vector evolves according 

to 

(3.11) v = Vm(# - /) = g ij ($ - f)iXj. 

Proof. Since v is unit normal vector we have v € T(M). Furthermore, differ- 
entiating 

(3.12) 0=(v,Xi) 
with respect to t, we deduce 

(3.13) (v,x i ) = -(v,x i ) = ($-j) i . □ 

3.4. Lemma (Evolution of the second fundamental form). The second fun- 
damental form evolves according to 

(3.14) h\ = (# - f)i + a{<£ - f)h!h{ + v($ - f)R aPlS v a x^x{g^ 
and 



(3.15) hij = (<P- f)ij - a(<P - fihfhkj + a($ - f)R a ^s 



a *i S 
V XiV'XA. 



Proof. We use the Ricci identities to interchange the covariant derivatives of v 
with respect to t and £ l 

(316) l^l^^li-n/^ 

- g kl (<P - f) ki x? + - faxfi - R a 0jS ^x]x s . 

For the second equality we used (|3 . 1 1(1 . On the other hand, in view of the 
Weingarten equation we obtain 

(3.17) £(i/f) = %(h$a%) = h\xt + h1± a k . 
Multiplying the resulting equation with g a px^ we conclude 

(3.18) h k g kj - a(0 - f)h k h k] = (* - f)ij + a($ - f)R aPlS v a x l3 l ^x s 3 

or equivalently (|3.14p . 

To derive (|3.15p , we differentiate 

(3.19) hij = h$g kj 

with respect to t and use (|3.8|) . □ 
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We emphasize that equation (|3.14p describes the evolution of the second 
fundamental form more meaningfully than (|3.15|) . since the mixed tensor is 
independent of the metric. 

3.5. Lemma (Evolution of (#— /)). The term (<i>— /) evolves according to 
the equation 



(# - /)' - ~ f)ij = rr./>f - /) + aj a v a {<f - /) 

+ <jSF^R a ^ a x^^x 5 J<P - /), 



(3.20) 
where 

(3.21) (^_/y = |(#-/) 

and 

(3.22) <J> = -^$(r). 

dr 

Proof. When we differentiate F with respect to t we consider F to depend on 
the mixed tensor h\ and conclude 

(3.23) (<P-fY = <PF;hl-f a i a ; 

The equation ([3~2"Uf then follows in view of ([3~T]) and (|3~T4"|) . □ 

3.6. Remark. The preceding conclusions, except Lemma 13.51 remain valid 
for flows which do not depend on the curvature, i.e., for flows 

x = —<j{—f)v = crfv. 

(3-24) 

x{0) = x , 

where / = f(x) is defined in an open set fl containing the initial spacelikc 
hypersurface Mq. In the preceding equations we only have to set = and 
/=/• 

The evolution equation for the mean curvature then looks like 

(3.25) H = -Af -a{\A\ 2 +R af3 v a ^}f, 

where the Laplacian is the Laplace operator on the hypersurface M(t). This is 
exactly the derivative of the mean curvature operator with respect to normal 
variations as we shall see in a moment. 

But first let us consider the following example. 

3.7. Example. Let (x a ) be a future directed Gaussian coordinate system 
in N, such that the metric can be expressed in the form 

(3.26) ds 2 = e 2 ^{a{dx Q ) 2 + a zj dx l dx j }. 

Denote by M(t) the coordinate slices {x° = t}, then M(t) can be looked at as 
the flow hypersurfaces of the flow 

(3.27) x = -a{-e^)v, 
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where we denote the geometric quantities of the slices by gij, v, hij, etc. 
Here x is the embedding 

(3.28) x = x{t,C) = {t,x i ). 

Notice that, if N is Riemannian, the coordinate system and the normal are 
always chosen such that v° > 0, while, if N is Lorentzian, we always pick the 
past directed normal. 

Hence the mean curvature of the slices evolves according to 

(3.29) H = -Ae* - a{\A\ 2 + R^v 13 }^ . 

We can now derive the linearization of the mean curvature operator of a 
spacelike hypersurface, compact or non-compact. 



3.8. Let Mq C N be a spacelike hypersurface of class C 4 . We first assume 
that Mo is compact; then there exists a tubular neighbourhood U and a cor- 
responding normal Gaussian coordinate system (x a ) of class C 3 such that 

is normal to Mo. 

Let us consider in U of Mo spacelike hypersurfaces M that can be writ- 
ten as graphs over Mo, M = graph u, in the corresponding normal Gaussian 
coordinate system. Then the mean curvature of M can be expressed as 

(3.30) H = {-Au + H-av- 2 u i u j h ij }v, 

where a — (is, v), and hence, choosing u = eip, ip £ C 2 (Mq), we deduce 

(3.31) de 

= -A v -o{\A\ 2 + R a0 u a ^)ip, 

in view of ([3~29]) . 

The right-hand side is the derivative of the mean curvature operator applied 

tO if. 

If Mo is non-compact, tubular neighbourhoods exist locally and the relation 
(|3.3ip will be valid for any ip £ C 2 (Mq) by using a partition of unity. 

The preceding linearization can be immediately generalized to a hypersurface 
Mo solving the equation 

(3-32) F, Mo = /, 

where / = f(x) is defined in a neighbourhood of Mo and F — F(hij) is curva- 
ture operator. 

3.9. Lemma. Let M be of class C m ' Q , m > 2, < a < 1, satisfy (|3732j) . 
Let U be a (local) tubular neighbourhood of Mo, then the linearization of the 
operator F — f expressed in the normal Gaussian coordinate system (x a ) cor- 
responding to U and evaluated at Mq has the form 

(3.33) - F ij Uij - n\F'*h l ,h,, + F ij R a p jS v a xfv' 1 x 5 j + f a v a }u. 
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where u is a function defined in Mq, and all geometric quantities are those of 
Mo; the derivatives are covariant derivatives with respect to the induced metric 
of Mq. The operator will be self-adjoint, if F lJ is divergence free. 

Proof. For simplicity assume that Mq is compact, and let u G C 2 (Mo) be fixed. 
Then the hypersurfaces 

(3.34) M e = graph(eu) 

stay in the tubular neighbourhood 11 for small e, |e| < eo, and their second 
fundamental forms (hij) can be expressed as 

(3.35) v~ x hij = —(eu)ij + hij, 

where hij is the second fundamental form of the coordinate slices {a; = const}. 
We are interested in 

(3-36) Te {F ~ f) ^ 

To differentiate F with respect to e it is best to consider the mixed form 
(h£) of the second fundamental form to derive 

(3 - 37) i {F ~ f) = F M ~§^ u = - Fi3u « + F & u - ^ 



where the equation is evaluated at e = and h\ is the derivative of h\ with 
respect to x°. 

The result then follows from the evolution equation (|3 . 14(1 for the flow ()3.27p , 
i.e., we have to replace (<£ — /) in (|3.14p by —1. □ 

4. Essential parabolic flow equations 

From (|3.14p on page [5] we deduce with the help of the Ricci identities a 
parabolic equation for the second fundamental form 

4.1. Lemma. The mixed tensor h\ satisfies the parabolic equation 
H-*F k % kl = 

o-$F M hrkh\hi - a<PFh„h rj + a(<P - f)h\h{ 

- f a px?49 kj + of a v a h\ + $F kl ' Ts h m h r J 
+ $FiFi + 2SF kl R af3jS x^xlx 5 r hrg rj 

- $F kl R a ^x a m x{x~lx\h™g T i - ^R^sx^x^xfh^ 
+ a$F kl R af3jS v a x l3 k v-<x s l hl - aSFR a ^ a xf^x 5 m g mj 

+ SF kl R a ^sAv a 4xl4x e m 9 mJ + v*x?xZx s m x' ig m 1}. 



(4.1) 
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Proof. We start with equation (|3.14p on page E] and shall evaluate the term 
(4-2) (*-f){\ 

since we are only working with covariant spatial derivatives in the subsequent 
proof, we may — and shall — consider the covariant form of the tensor 

(4-3) (*-/)«• 
First we have 

(4.4) # 4 = m = SF kl h k i,i 
and 

(4.5) <% = $F kl hki;ij + SF kl h k i, t F rs h rs . 3 + •i>F u -h u ,,h r>:j . 

Next, we want to replace hki-ij by h^-ki ■ Differentiating the Codazzi equation 

(4.6) hki-,i = hik-i + R a fi~ 1 bv a x li k x{x\, 
where we also used the symmetry of hik, yields 

hfcl-ij — hik-ij + R a ^jS;e^ x k x l x i x j 

+ RaPiSiVj Vk X l X i + V X kj X l x i + v x k x lj x i + v x k x l x ij}- 

To replace hki-ij by hij-^i we use the Ricci identities 

(4.8) hik-ij = hik-ji + h a kR a uj + h ai R a kl j 
and differentiate once again the Codazzi equation 

(4.9) hik-j = hij-k + Rupys^x^xlxj. 
To replace we use the chain rule 

ft fa%i j 



(4-10) 



fij — fa0 x i Xj + faX^j- 



Then, because of the GauB equation, Gaussian formula, and Weingarten 
equation, the symmetry properties of the Riemann curvature tensor and the 
assumed homogeneity of F, i.e., 

(4.11) F = F kl h ku 

we deduce ||HT) from (|3TT4j) on page [6] after reverting to the mixed representa- 
tion. □ 

4.2. Remark. If we had assumed F to be homogeneous of degree do instead 
of 1, then we would have to replace the explicit term F — occurring twice in the 
preceding lemma — by dyF. 

If the ambient semi-Ricmannian manifold is a space of constant curvature, 
then the evolution equation of the second fundamental form simplifies consid- 
erably, as can be easily verified. 
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4.3. Lemma. Let N be a space of constant curvature Kn , then the second 
fundamental form of the curvature flow (|3.7|) on page [5] satisfies the parabolic 
equation 

hi - <PF kl hi M = a$F kl h rk h\h\ - a$Fh„h r J + o(<P - ])h\h? k 
- Upxfx^ + u] a v a h\ + &F kl > rs h m h rs i 
+ <PFiF J 

+ K N {($ - f)5j + SF6{ - SF kl g kl hj}. 

Let us now assume that the open set fi C N containing the flow hyper- 
surfaces can be covered by a Gaussian coordinate system (x a ), i.e., Q can be 
topologically viewed as a subset of / x iSo, where So is a compact Riemannian 
manifold and / an interval. We assume furthermore, that the flow hypersur- 
faces can be written as graphs over Sq 

(4.13) M(t) = { x° = u(x l ) : x = (x 1 ) G S Q }; 

we use the symbol x ambiguously by denoting points p = (x a ) S N as well as 
points p = (x l ) G So simply by x, however, we are careful to avoid confusions. 

Suppose that the flow hypcrsurfaces are given by an embedding x — x(t, £), 
where £ = are local coordinates of a compact manifold Mo, which then has 
to be homeomorphic to So, then 

x°=u(t, 0=u{t,x(t,0), 

(4.14) 

x l =x\t,i). 
The induced metric can be expressed as 

(4.15) gij = (xi,Xj) = auiU 3 +o- kl x k x l J , 
where 

(4.16) u t = u k x k , 
i.e., 

(4.17) = {cru k ui + o k {\x k x l j, 

hence the (time dependent) Jacobian (x k ) is invertible, and the (£*) can also 
be viewed as coordinates for So- 

Looking at the component a = of the flow equation (|3.7[) on page [5] we 
obtain a scalar flow equation 

(4.18) ii= -e^v^i-P- f), 

which is the same in the Lorentzian as well as in the Riemannian case, where 

(4.19) v 2 = 1 - aa ij UiUj, 
and where 

(4.20) \Du\ 2 = a l3 u lUj 
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is of course a scalar, i.e., we obtain the same expression regardless, if we use 
the coordinates x % or £ 4 . 

The time derivative in (|4.18p is a total time derivative, if we consider u to 
depend on u = u(t, x(t, £,)). For the partial time derivative we obtain 



where u l = a iJ Uj. 

Controlling the C 1 -norm of the graphs M(t) is tantamount to controlling v, 
if N is Riemannian, and v = v , if N is Lorentzian. The evolution equations 
satisfied by these quantities are also very important, since they are used for 
the a priori estimates of the second fundamental form. 

Let us start with the Lorentzian case. 

4.4. Lemma (Evolution of v). Consider the flow (|3.7p in a Lorentzian space 
N such that the spacelike flow hyper surfaces can be written as graphs over Sq. 
Then, v satisfies the evolution equation 



where r\ is the covariant vector field {rj a ) = e^(— 1, 0, . . . , 0). 

Proof. We have v = (r],v). Let be local coordinates for M(t). Differenti- 
ating v covariantly we deduce 




(4.23) 



$F^h lk h k v + [(# - /) - $F]r, af3 v a v 
2£F i Jrfx?x p k r ]a0 - £F i i Va frX?x]v a 
SF i ^R aM& v a xP i xlx 5 j r le xtg kl 



(4.24) 



Vi = r\ a f}X l3 i v a + r\ a v\ 



,9 
i > 



Vij = ■q al 3 1 x'lx1v a + 7y Q/ 3Ccfi/ Q 
25) 3 v 

+ Vaptf v <* + r\ afi x p -v°; + rj a v\ 

The time derivative of v can be expressed as 



v = rj a p± l3 v' x + rj a v 



(4.26) 



r, a pis a ^(<P-f) + (<P-f) k x k <r ] , 



la 



r^z/V^ - /) + SF k xZ Va - f x^xtg ik v, 



where we have used p. lip on page [6] 
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Substituting (|4.25p and ()4.26|1 in (|4.23|1 . and simplifying the resulting equa- 
tion with the help of the Weingarten and Codazzi equations, we arrive at the 
desired conclusion. □ 

In the Riemannian case we consider a normal Gaussian coordinate system 
(x a ), for otherwise we won't obtain a priori estimates for v, at least not without 
additional strong assumptions. We also refer to x° = r as the radial distance 
function. 

4.5. Lemma (Evolution of v). Consider the flow (|3.7p in a normal Gaussian 
coordinate system where the M(t) can be written as graphs of a function u(t) 
over some compact Riemannian manifold Sq . Then the quantity 



(4.27) v = y/l + \Du\ 2 = (r a v a ) 

satisfies the evolution equation 



(4.28) 



+ r a/3 u a vP[($ - f) - <PF}v 2 + 2$F ij h k i r otf} x%x?v 1 

+ QF^R^^x^xlr^g^v 2 

+ $F*r a ( h v a xfx]i? + ~f a x a m g mk rpxlv 2 . 



Proof. Similar to the proof of the previous lemma. □ 

The previous problems can be generalized to the case when the right-hand 
side / is not only defined in TV or in Q but in the tangent bundle T(N) resp. 
T(f2). Notice that the tangent bundle is a manifold of dimension 2(n + 1), i.e., 
in a local trivialization of T(N) f can be expressed in the form 

(4.29) f = f{x,u) 

with x e N and v € T X (N), cf. 17, Note 12.2.14]. Thus, the case / = f{x) is 
included in this general set up. The symbol v indicates that in an equation 

(4-30) F, M =/(z,i/) 

we want / to be evaluated at (x,u), where x G M and v is the normal of M in 
x. 

The Minkowski problem or Minkowski type problems are also covered by the 
present setting, though the Minkowski problem has the additional property that 
the problem is transformed via the Gaufl map to a different semi-Riemannian 
manifold as a dual problem and solved there. Minkowski type problems have 
been treated in [S], [53], [H] and [3T]- 

4.6. Remark. The equation (|4.30p will be solved by the same methods as 
in the special case when / = f(x), i.e., we consider the same curvature flow, 
the evolution equation (|3.7|) on pageO as before. 
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The resulting evolution equations are identical with the natural exception, 
that, when / or / has to be differentiated, the additional argument has to be 
considered, e.g., 

(4.31) h = fax? + M = hx a + f„ e 4 h\ 
and 

(4.32) / = f a x a + fait = -a($ - j)f a u a + f„ e g ij ($ - /)<af. 
The most important evolution equations are explicitly stated below. 

Let us first state the evolution equation for {<P — /) . 

4.7. Lemma (Evolution of {<&— /)). The term (<P — f) evolves according to 
the equation 

(<P - /)' - &F^(<P - /) y = a.N^h, l .h l -:<l> - /) 

(4.33) + *f a v a ($ - f) - /„«<(# - f) j9 ij 

+ a^R a&1& v a x^x\(^ - /), 

where 

(4-34) (*_/)' = fL(*_/) 



(4.35) $=4-$(r). 

dr 

Here is the evolution equation for the second fundamental form. 
4.8. Lemma. The mixed tensor h\ satisfies the parabolic equation 

K-i>F k % kl 

= G^F kl h rk h\h\ - a<PFh rl h ri + a($ - f)h k h{ 

- Upxfx^ + oj a v a h\ - f a ^(xTx k h k > +x?xlh k g lj ) 

- f vav ex?xlh k h l i - L,x k h k t g l > + af v „u a h k hi 
(4-36) + $ F ki,rs hki . Ks i + 2<PF kl R a p lS x^xfxlx 5 r hrg r3 

- $F kl R a p l8 x a m xlx~lx\h™g r i - ^R^sx^x^xfh^ 
+ o^F kl R aM&V a xlv^x\h\ - o-&FR af37 sv a xf^x s m g mj 

+ <j(<P - f)R a p jS v a x?^x 5 m g m i + &FiF j 

+ <PF kl R a ^sA^x^xtx^ + v a x?xlx 5 m xtg mj }. 
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The proof is identical to that of Lemma 14. 1| we only have to keep in mind 
that / now also depends on the normal. 

If we had assumed F to be homogeneous of degree do instead of 1, then, we 
would have to replace the explicit term F — occurring twice in the preceding 
lemma — by d^F . 



4.9. Lemma (Evolution of v). Consider the flow (|3.7|) in a Lorentzian space 
N such that the spacelike flow hypersurfaces can be written as graphs over Sq . 
Then, v satisfies the evolution equation 



v - SF^Vij = - '/'/ ''//,i//;r + [(# - /) - SF] Va0 u a ^ 

(4.37) 



2<PF^h k 3 xfxli la p - &FVr) a p 7 x?x]v a 



- &F i iR af3l sv a xfxlx 5 j r }e x* l g kl 
~ fp^x^n a g ik - f u0 x%h ik x?ri a , 

where n is the covariant vector field (r] a ) = e^(— 1, 0, . . . , 0). 



The proof is identical to the proof of Lemma 14.41 
In the Riemannian case we have: 



4.10. Lemma (Evolution of v). Consider the flow Q3.7P in a normal Gauss- 
ian coordinate system (x a ), where the M(t) can be written as graphs of a func- 
tion u(t) over some compact Riemannian manifold Sq. Then the quantity 



(4.38) v = ^l + \Du\ 2 = {r a v a )~ 1 
satisfies the evolution equation 

v - QF^Vij = - -/>/ ''//, a //"Y - 2v- 1 $F lj v l v j 
+ [(# - /) - SF]r a pv a ^v 2 

(4.39) + 2$F^h k x?x k r a pv 2 + 4>F ij r af } 7 x?x]v a v 2 

+ &F i iR a p yS v a xfx'lx s j r e x'ig kl v 2 

+ f fi x lX tr a g lk v 2 + f^x^h tk xfr a v 2 , 
where r — x° and (r a ) — (1, 0, . . . , 0). 



5. Stability of the limit hypersurfaces 

5.1. Definition. Let N be semi-Riemannian, F a curvature operator, and 
M C N a compact, spacelike hypersurface, such that M is admissible and 
F lJ , evaluated at (hij,gij), the second fundamental form and metric of M, is 
divergence free, then M is said to be a stable solution of the equation 



(5.1) 



Ft = f 

M J 1 
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where / = f(x) is defined in a neighbourhood of M, if the first eigenvalue Ai 
of the linearization, which is the operator in (I3.33|) on pageJH is non- negative, 
or equivalently, if the quadratic form 

(5.2) f F^muj -a I {^hfhkj + F ij R a ^ a x^"'x 5 j + f a v a }u 2 
Jm Jm 

is non-negative for all u £ C 2 (M). 

It is well-known that the corresponding eigenspace is then onedimensional 
and spanned by a strictly positive eigenfunction 77 

(5.3) - F^rhj - a{F^h^h kj + F 1 ^ R aM& v a x^ x] + f a v a }n = X lV . 

Notice that F^ is supposed to be divergence free, which will be the case, if 
F = Hk, 1 < k < n, and the ambient space has constant curvature, as we 
shall prove at the end of this section. If k = 1 , then F^ = g % i and N can be 
arbitrary while in case k — 2, we have 

(5.4) /•" ////' I,' 1 . 

hence N Einstein will suffice. 

To simplify the formulation of the assumptions let us define: 

5.2. Definition. A curvature function F is said to be of class (D), if for 
every admissible hypersurface M the tensor F lJ , evaluated at M, is divergence 
free. 

We shall prove in this section that the limit hypersurface of a converging 
curvature flow will be a stable stationary solution, if the initial flow velocity 
has a weak sign. 

5.3. Theorem. Suppose that the curvature flow (|3.7p on page [S] exists for 
all time, and that the leaves M(t) converge in C 4 to a hypersurface M , where 
the curvature function F is supposed to be of class (D) . Then M is a stable 
solution of the equation 

(5-5) F lu = f 

provided the velocity of the flow has a weak sign 

(5.6) <?-/>0 V <2>-/<0 

at t = and M (0) is not already a solution of (|5.5p . 

Proof. Convergence of a subsequence of the M(t) would actually suffice for 
the proof, however, the assumption (|5.6p immediately implies that the flow 
converges, if a subsequence converges and a priori estimates in C 4 ' Q are valid. 

The starting point is the evolution equation (|3.20[) on page [7] from which we 
deduce in view of the parabolic maximum principle that <P — f has a weak sign 
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during the evolution, cf. |18[ Proposition 2.7.1], i.e., if we assume without loss 
of generality that at t = 

(5.7) # - / > 0, 

then this inequality will be valid for all t. Moreover, there holds 

(5.8) f (<2>-/)>0 V0<t<oo 

J M(t) 

if this relation is valid for t — 0, as we shall prove in the lemma below. 
On the other hand, the assumption 



(5-9) / (*-/)> 

Jm(o) 

is a natural assumption, for otherwise the initial hypersurface would already 
be a stationary solution which of course may not be stable. 

Notice also that apart from the factor <f> the equation (|3.20[) looks like the 
parabolic version of the linearization of (F — /). If the technical function 
<P = <P(r) is not the trivial one <P(r) = r, then we always assume that / > 
and that this is also valid for the limit hypersurface M. Only in case <P(r) = r 
and F — H, we allow / to be arbitrary. 

Thus, our assumptions imply that in any case 

(5.10) <i>>e o >0 VieIR + . 

Furthermore, we derive from (|3.20p that not only the elliptic part converges 
to but also 

(5.11) (* - /)' = <PF + a<P(f)f a ^ a (<P - /), 



(5.12) limF = 0. 

Suppose now that M is not stable, then the first eigenvalue Ai is negative 
and there exists a strictly positive eigenfunction r/ solving the equation (|5-3[) 
evaluated at M. Let U be a tubular neighbourhood of M with a corresponding 
future directed normal Gaussian coordinate system (x a ) and extend 77 to U by 
setting 

(5.13) v{x a ,x) = r)(x), 

where, by a slight abuse of notation, we also denote (x l ) by x. Thus there holds 

(5.14) r) a v a = 

in M, and choosing U sufficiently small, we may assume 

(5.15) \v a v a \<V 

for all hypersurfaces M(t) C U. 
Now consider the term 

(5.16) / 

J M(t) 
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for large t, which converges to 0. Since it is positive, in view of (|5.8|) . there 
must exist a sequence of t, not explicitly labelled, tending to infinity such that 



at JM(t) 

(5.17) = / -$- 2 $F{$ - f)n + / <i>- 1 (<£-/)'r/ 

JM(t) JM(i) 

-a I &~ l ri*v a ($ ~ ff -tr f f) 2 Hr h 

JM(t) J M(t) 

where we used the relation ()3.8|) on page [5] to derive the last integral. 

The rest of the proof is straight-forward. Multiply the equation (|3 . 20(1 by 
l P^ 1 r] and integrate over M(t) for those values of t satisfying the preceding 
inequality to deduce 



(5.18) 



JM(t) JM(t) 



a [ {F^h\h k , + F^R^v^x^x] + $- 1 $(J)f a v a }ri($ - /), 

JM(t) 

and conclude further that the right-hand side can be estimated from above by 

(5.19) ^/ T,($-j) 

1 JM(t) 

for large t, while the left-hand side can be estimated from below by 



(5.20) -e{t) (<P-f)v 

JM(t) 

such that 

(5.21) e(t)>0 A lime(t) = 

in view of (|5.17|) , where we used (|5.12p . (|5.15[) as well as 

(5.22) Iim(* - /) = 0; 

a contradiction because of (15.81). □ 



5.4. Lemma. Let M(t) be a solution of the curvature flow (|3.7p on page [5] 
defined on a maximal time interval [0, T*), < T* < oo, and suppose that 
<P — f has a weak sign at t = 0, e.g., 

(5.23) (# - /) > 
and suppose furthermore that 

(5.24) f (<2> - /) > 0, 

Jm(o) 

then 

(5.25) ( (<P-/)>0 V0<i<T*. 

JM(t) 
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Proof. Let Mq be an abstract compact Riemannian manifold that is being 
isometrically embedded in N with image M(0). Let be a generic coordinate 
system for Mq and abbreviate — /) by u. The evolution equation ()3.20p can 
then be looked at as a linear parabolic equation for u — u(t,£) on Mq with 
time dependent coefficients and time dependent Riemannian metric <7ij (£,£)• 

By assumption u doesn't vanish identically at t — 0, i.e., there exists a ball 
Bp = -Bp (Co) such that 

(5.26) «(0,f)>° VeeB p (Co). 
Let C be the cylinder 

(5.27) C = [0,T*)xB p 
and assume that there exists a first to > such that 

(5.28) infu(fo,-) = = u(*o,fi). 

b p 

We shall show that this is not possible: If £t g then this contradicts the 
strong parabolic maximum principle, cf. [ISJ Lemma 2.7.1], and if £i S 9i? p , 
then we deduce from [lSl Lemma 2.7.4] (a parabolic version of the Hopf Lemma) 

Fin 

(5.29) £^ , 6 ) <0 , 

where is the exterior normal of the ball B p in £i, contradicting the fact that 
the gradient of u(toj ■) vanishes in £i because it is a minimum point; notice that 
we already know u > in [0, T*) x M . □ 

For some curvature operators one can prove a priori estimates for the second 
fundamental form only for stationary solutions and not for the leaves of a 
corresponding curvature flow. In order to use a curvature flow to obtain a 
stationary solution one uses „e-regularization", i.e., instead of the curvature 
function F one considers 

(5.30) F(h lJ )=F(h lJ +eHg lJ ) 

for e > 0, and starts a curvature flow with F and fixed e > 0. 

A priori estimates for the regularized flow are usually fairly easily derived, 
since 

(5.31) /' '•' /•"' • </' 

but of course the estimates depend on e. Having uniform estimates one can 
deduce that the flow — or at least a subsequence — converges to a limit hyper- 
surfaces M e satisfying 

(5-32) F lMc = f. 

Then, if uniform C 4, "-estimates for the M e can be derived, a subsequence will 
converge to a solution M of 

(5.33) F> M = /, 
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cf. [13] , where this method has been used to find hypersurfaccs of prescribed 
scalar curvature in Lorentzian manifolds, see also Theorem 16.91 on page 1361 

We shall now show that the solutions M obtained by this approach are all 
stable, if F is of class (D) and the initial velocities of the regularized flows have 
a weak sign. Notice that the curvature functions F are in general not of class 
(D). 

5.5. Theorem. Let F be of class (D), then any solution M of 
(5-34) F lM = f 

obtained by a regularized curvature flow as described above is stable, provided 
the initial velocity of the regularized flow has a weak sign, i.e., it satisfies 

(5.35) <5 - / > V <2> - / < 

at t = and the flow hypersurfaccs converge to the stationary solution in C 4 . 

Proof. Let M e be the limit hypersurfaces of the regularized flow for e > 0, and 
assume that the M e satisfy uniform C '"-estimates such that a subsequence, 
not relabelled, converges in C 4 to a compact spacelike hypersurface M solving 
the equation 

(5.36) F {M = f. 

Assume that M is not stable so that the first eigenvalue of the linearization 
is negative and there exists a strictly positive eigenfunction n satisfying (|5.3|) . 
Extend 77 in a small tubular neighbourhood U of M such that (|5.15j) is valid 
for all M e , if e is small, e < eo- 

For those e we then deduce 

- vi&hfhv + Pv R^u^u^x] + f a v a }r) < ^77, 

where the inequality is evaluated at M e and where we used the convergence in 
C 4 . 

Now, fix e, e < eo, then the preceding inequality is also valid for the flow 
hypersurfaces M(t) converging to M e , if t is large, and the same arguments as 
at the end of the proof of Theorem 15.31 lead to a contradiction. Hence, M has 
to be a stable solution. □ 

Knowing that a solution is stable often allows to deduce further geometric 
properties of the underlying hypersurface like that it is either strictly stable 
or totally geodesic especially if the curvature function is the mean curvature, 
cf. e.g., [25], where the stability property has been extensively used to deduce 
geometric properties. 

We want to prove that a neighbourhood of stable solutions can be foliated 
by a family of hypersurfaces satisfying the equation modulo a constant. 
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5.6. Theorem. Let M (Z N be compact, spacelike, orientable and a stable 
solution of 

(5-38) G| M = (F-/)| M =0, 

where F is of class (D) and M as well as F, f are of class C m - a , 2 < m < oo, 
< a < 1, then a neighbourhood of M can be foliated by a family 

(5.39) A = {M e : \e\ < e } 
of spacelike C m ' a -hyper surf aces satisfying 

(5.40) G| Me =r(e), 

where r is a real function of class C m ' a . The M e can be written as graphs over 
M in a tubular neighbourhood of M 

(5.41) M f = { (u(e, x),x): x G M} 

such that u is of class C m ' a in both variables and there holds 

(5.42) u > 0. 

Proof, (i) Let us assume that M is strictly stable. Consider a tubular neigh- 
bourhood of M with corresponding normal Gaussian coordinates (x a ) such 
that M — {x° = 0}. The nonlinear operator G can then be viewed as an 
elliptic operator 

(5.43) G : B p (0) C C m ' a (M) -» C m -^ a {M) 

where p is so small that all corresponding graphs are admissible. 

In a smaller ball DG is a topological isomorphism, since M is strictly stable, 
and hence G is a diffeomorphism in a neighbourhood of the origin, and there 
exist smooth unique solutions 

(5.44) M e = {u(e,x):xeM} \e\ < e 
of the equations 

(5.45) G lMc = e 

such that u e C m ' a ((-e ,e ) x M). 
Differentiating with respect to e yields 

(5.46) DGii = 1. 

Let us consider this equation for e = 0, i.e., on M, and define 

(5.47) f] = min(w, 0). 
Then we deduce 

(5.48) 0< [ (DGr], rj) = [ r, < 0, 

Jm Jm 

and hence there holds 

(5.49) u > 0, 
because of the strict stability of M. 
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Applying then the maximum principle to (15.46|) . we deduce further 

(5.50) infu>0, 

M 

hence the hypersurfaces form a foliation if eo is chosen small enough such that 

(5.51) infw(e,-)>0 V|e|<e . 

(ii) Assume now that M is not strictly stable. After introducing coordinates 
corresponding to a tubular neighbourhood hi of M as in part (i) any function 
u € C' m,a (M) with |u| m . a small enough defines an admissible hypersurface 

(5.52) M(u) = graphic U 
such that G\ M . u) can be expressed as 

(5.53) G| M(u) = G(u). 
Let 

(5.54) A = DG(0), 
then A is self-adjoint, monotone 

(5.55) (Au,u)>0 VueH l < 2 (M) 

and the smallest eigenvalue of A is equal to zero, the corresponding eigenspace 
spanned by a strictly positive eigenfunction r\. 

Similarly as in [3J p. 621] we consider the operator 

(5.56) *(u,r) = (G(u) -r,ip(u)) 

defined in B p (0) x R, B p (Q) C C m > a (M) for small p > 0, where ip is a linear 
functional 

(5.57) L p{ u ) = i] u 

Jm 

*f? is of class C m ' a and maps 

(5.58) * : B p (0) xR^ C m ^ a (M) x M, 
such that 



(5.59) ZM 



£>G -1\ 

evaluated at (0, 0) is bijective as one easily checks. Indeed let (u, e) satisfy 

(5.60) D*(u,e) = (0,0), 
then 

(5.61) Au = DGu = € A / r)u = 0, 



hence 



(5.62) e / 7j = (Au, r?) = (it, A77) = 
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and we conclude e = as well as u = 0. 

To prove the surjectivity, let (w, S) £ C m ~ 2,a (M) xl be arbitrary. Choosing 



(5.63) 



Jm W7 1 



we deduce 

(5.64) [ (e + w)r/ = 0, 

Jm 

hence there exists u £ C m,a (M) solving 

(5.65) Au = e + w 
and 

(5.66) u = u + Xrj 
with 

(5.67) \ = 5- I rju 

JM 

then satisfies 

(5.68) / rju = e, 
i.e., 

(5.69) DV(u,e) = (w,5). 

Applying the inverse function theorem we conclude that there exists eo > 
and functions (u(e, x), r(e)) of class C m,a in both variables such that 

(5.70) G(u(e))=T(e) A / rju{e) = e V|e|<e ; 

Jm 

r(e) is constant for fixed e. 
The hypersurfaces 

(5.71) A = {M e = M(u(e)): |e| < e } 
will form a foliation, if we can show that 

(5.72) u ± 0. 

Differentiating the equations in (|5.70|) with respect to e and evaluating the 
result at e = yields 

(5.73) Au(Q) = f(0) A / 776(0) = 1 

Jm 

and we deduce further 

(5.74) f (0) f r) = (Au(0),r)) = (u(0),A V ) = 

Jm 

and thus 

(5.75) f(0) = A u(0) = 77 > 0, 
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if r\ is normalized such that (77,77) = 1, i.e., we have u{e) > 0, if eo is chosen 
small enough. □ 

5.7. Remark. Let M be a stable solution of 

(5.76) G\ M = 

as in the preceding theorem, but not strictly stable and let M t be a foliation 
of a neighbourhood of M such that 

(5.77) G\ Mc =r(e) V| £ |<e . 

If M is the limit hypersurface of a curvature flow as in Theorem [531 then 

(5.78) r(e)>0 V0 < e < e , 

if the flow hypersurfaces M(t) converge to M from above, which is tantamount 
to 

(5.79) <P(F) - / > 0, 
or we have 

(5.80) r(e)<0 V - e < e < 0, 
if 

(5.81) 0(F) - f < 0, 

in which case the flow hypersurfaces converge to M from below. 

The direction „ above" is defined by the region the normal av of M points 

to. 

Proof. Let us assume that the flow hypersurfaces satisfy (15.79P and fix < e < 
eo. We may also suppose that the initial hypersurface M(0) doesn't intersect 
the tubular neighbourhood of M which is being foliated by M e . Now, fix 
< e < eo, then there must be a first t > such that M(t) touches M e from 
above which yields, in view of the maximum principle, 

(5.82) G, Mc = r(e) > 0, 

since r(e) < would imply r(e) = and M e = M(t), cf. [H Theorem 2.7.9], 
i.e., M(t) would be a stationary solution, which is impossible as we have proved 
in Lemma [BTH □ 

Finally, let us show that the symmetric polynomials Hk, 1 < k < n, are of 
class (D), if the ambient space has constant curvature. 

5.8. Lemma. Let N be a semi-Riemannian space of constant curvature, 
then the symmetric polynomials F = Hk. 1 < k < n, are of class (D). In case 
k = 2 it suffices to assume N Einstein. 
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Proof. We shall prove the result by induction on k. First we note that the 
cones of definition Jfc C M™ of the Hj- form an ordered chain 

(5.83) r k C A-i VKKn, 

cf. [7], so that a hypersurface admissible for Hk is also admissible for H^-i- 
For fc = 1 we have 

(5.84) F ij = g 13 

and the result is obviously valid for arbitrary N. 

Thus let us assume that the result is already proved for 1 < k < n. Set 
F = Hk+i, F = Hk and let M be an admissible hypersurface for F with 
principal curvatures Kj. 

From the definition of the H^s we immediately deduce 

dF OF 

(5 - 85) F= W- +Ki W- 
for fixed i, no summation over i, or equivalently, 

(5.86) Fg^ =F l i +F jm V m , 

notice that the last term is a symmetric tensor, since for any symmetric cur- 
vature function F F' J and hij commute, cf. [TH1 Lemma 2.1.9]. Thus there 
holds 

(5.87) F l] = Fg 11 - F^ m h l m 
and we deduce, using the induction hypothesis, 

F.7 = F l - F im hL.j = F l - P m h ml l 

(5.88) J A A ' J J ' 

= F l - F l = 0, 

where we applied the Codazzi equations at one point. 
IfF = if 2 , then 

(5.89) F lj = Hg lj - h ij 

and the assumption N Einstein suffices to conclude that is divergence 
free. □ 



6. Existence results 

From now on we shall assume that ambient manifold N is Lorentzian, or 
more precisely, that it is smooth, globally hyperbolic with a compact, connected 
Cauchy hypersurface. Then there exists a smooth future oriented time function 
x° such that the metric in N can be expressed in Gaussian coordinates (x a ) as 

(6.1) ds 2 = e 2,p {-(dx ) 2 + a ij dx i dx j }, 

where x° is the time function and the (x % ) are local coordinates for 

(6.2) So = {x° = 0}. 
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Sq is then also a compact, connected Cauchy hypersurface. For a proof of the 
splitting result see [H Theorem 1.1], and for the fact that all Cauchy hyper- 
surfaces are diffeomorphic and hence Sq is also compact and connected, see [3J 
Lemma 2.2]. 

One advantage of working in globally hyperbolic spacetimes with a compact 
Cauchy hypersurface is that all compact, connected spacelike C m -hypersurfaces 
M can be written as graphs over Sq. 

6.1. Lemma. Let N be as above and M C N a connected, spacelike hyper- 
surface of class C rn , 1 < m, then M can be written as a graph over So 

(6.3) M = graphu| 5o 

with u G C m (<S ). 

We proved this lemma under the additional hypothesis that M is achronal, 
|10l Proposition 2.5], however, this assumption is unnecessary as has been 
shown in [551 Theorem 1.1]. 

We are looking at the curvature flow (|3.7p on page [5] and want to prove that 
it converges to a stationary solution hypersurface, if certain assumptions are 
satisfied. 

The existence proof consists of four steps: 

(i) Existence on a maximal time interval [0, T*). 

(ii) Proof that the flow stays in a compact subset. 

(iii) Uniform a priori estimates in an appropriate function space, e.g., C 4 ' Q (So) 
or C° c (So), which, together with (ii), would imply T* — oo. 

(iv) Conclusion that the flow — or at least a subsequence of the flow hypersur- 
faces — converges if t tends to infinity. 

The existence on a maximal time interval is always guaranteed, if the data 
are sufficiently regular, since the problem is parabolic. If the flow hypersurfaces 
can be written as graphs in a Gaussian coordinate system, as will always be the 
case in a globally hyperbolic spacetime with a compact Cauchy hypersurface 
in view of Lemma 16. 11 the conditions are better than in the general case: 

6.2. Theorem. Let 4 < m G N and < a < 1, and assume the semi- 
Riemannian space N to be of class C m+2 ' a . Let the strictly monotone curvature 
function F, the functions f and & be of class C m,Q and let M G C m+2 ' a be 
an admissible compact, spacelike, connected, orientabl^ hypersurface. Then the 
curvature flow (|3.7[) on page [5] with initial hypersurface Mq exists in a maximal 
time interval [0, T*), < T* < oo, where in case that the flow hypersurfaces 
cannot be expressed as graphs they are supposed to be smooth, i.e, the conditions 
should be valid for arbitrary 4 < m G N in this case. 



Recall that oriented simply means there exists a continuous normal, which will always 
be the case in a globally hyperbolic spacetime. 
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A proof can be found in [HI Theorem 2.5.19, Lemma 2.6.1]. 

The second step, that the flow stays in a compact set, can only be achieved 
by barrier assumptions, cf. Definition 12.11 Thus, let Q c N be open and 
precompact such that dfl has exactly two components 

(6.4) dn = Mi U M 2 

where Mi is a lower barrier for the pair (F, /) and M 2 an upper barrier. More- 
over, Mi has to lie in the past of M2 

(6.5) M 1 C/~(M 2 ), 

cf. [HI Remark 2.7.8]. 

Then the flow hypersurfaces will always stay inside fi, if the initial hyper- 
surface Ma satisfies Mo C £2, [HI Theorem 2.7.9]. This result is also valid if 
Mo coincides with one the barriers, since then the velocity — /) has a weak 
sign and the flow moves into Q for small t, if it moves at all, and the arguments 
of the proof are applicable. 

In Lorcntzian manifolds the existence of barriers is associated with the pres- 
ence of past and future singularities. In globally hyperbolic spacetimes, when 
N is topologically a product 

(6.6) N = IxS , 

where / = (a, b), singularities can only occur, when the endpoints of the interval 
are approached. A singularity, if one exists, is called a crushing singularity, if 
the sectional curvatures become unbounded, i.e., 

(6.7) R a p lS R al3lS -► 00 

and such a singularity should provide a future resp. past barrier for the mean 
curvature function H . 



6.3. Definition. Let TV be a globally hyperbolic spacetime with compact 
Cauchy hypersurface Sq so that N can be written as a topological product 
N = I x So and its metric expressed as 

(6.8) ds 2 = e 2,p (-(dx ) 2 + a l3 {x° ,x)dx l dx 3 ). 

Here, x° is a globally defined future directed time function and (x l ) are lo- 
cal coordinates for Sq. N is said to have a future resp. past mean curvature 
barrier, if there are sequences M^ resp. M^T of closed, spacelike, admissible 
hypersurfaces such that 

(6.9) lim H\ =00 resp. lim H 1 = —00 
and 

(6.10) lim sup inf x° > x°(p) VpE N 
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resp. 

(6.11) liminf supx < x°(p) M p € iV, 

If one stipulates that the principal curvatures of the M£ resp. tend to 
plus resp. minus infinity, then these hypersurfaces could also serve as barriers 
for other curvature functions. The past barriers would most certainly be non- 
admissible for any curvature function except H . 

6.4. Remark. Notice that the assumptions (|6.9[) alone already implies (|6. 10[) 
resp. (|6.1ip . if either 

(6.12) limsup inf x° > a 
resp. 

(6.13) liminf sup x° < b 
where (a, b) = x°(N), or, if 

(6.14) R a 0V a ^>-A V{v,u) = -1. 
where A > 0. 

Proof. It suffices to prove that the relation (|6.10p is automatically satisfied 
under the assumptions (|6.12|) or ()6.14j) by switching the light cone and replacing 
x° by —x° in case of the past barrier. 
Fix k, and let 

(6.15) Tfe = infx°, 
then the coordinate slice 

(6.16) M Tk = {x° = r fe } 

touches Mfe from below in a point pk G Mk where = x°(j>k) and the maxi- 
mum principle yields that in that point 

(6-17) H lMTk >H ]Mk , 

hence, if k tends to infinity the points (j>k) cannot stay in a compact subset, 
i.e., 

(6.18) limsupxVfe) -> b 
or 

(6.19) limsupx°(pfc) — > a. 

We shall show that only (|6.18|) can be valid. The relation (|6.19[) evidently 
contradicts (|6.12|) . 
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In case the assumption (|6.14j) is valid, we consider a fixed coordinate slice 
Mo = {a; = const}, then all hypersurfaces Mk satisfying 

(6.20) H\ M <in£H A vnl < inf H 

have to lie in the future of Mo, cf. [T8l Lemma 4.7.1], hence the result. □ 

A future mean curvature barrier certainly represents a singularity, at least 
if N satisfies the condition 

(6.21) R afj v a vP >-A V(v,v) = -1 

where A > 0, because of the future timelike incompleteness, which is proved in 
P], and is a generalization of Hawking's earlier result for A = 0, [24] . But these 
singularities need not be crushing, cf. [HI Section 2] for a counterexample. 

The uniform a priori estimates for the flow hypersurfaces are the hardest 
part in any existence proof. When the flow hypersurfaces can be written as 
graphs it suffices to prove C 1 and C 2 estimates, namely, the induced metric 

(6.22) g lJ (t,Z) = (x l ,x ] ) 

where x — x(t,£) is a local embedding of the flow, should stay uniformly 
positive definite, i.e., there should exist positive constants c,, 1 < i < 2, such 
that 

(6.23) c igi j(0,0 < 9ij(t,0 < c 2 5„(0,0, 
or cquivalcntly, that the quantity 

(6.24) v = {r),v), 

where v is the past directed normal of M(t) and r\ the vector field 

(6.25) ry = (7 7Q )=e^(-l,0,...,0) ) 

is uniformly bounded, which is achieved with the help of the parabolic equation 
(|4.37[) on page[l5j if it is possible at all. 

However, in some special situations C^-estimates are automatically satisfied, 
cf. Theorem 16. Ill at the end of this section. 

For the C 2 -estimates the principle curvatures k.; of the flow hypersurfaces 
have to stay in a compact set in the cone of definition F of F, e.g., if F is the 
Gaussian curvature, then r = _T + and one has to prove that there are positive 
constants fcj, i = 1, 2 such that 

(6.26) ki<Ki<k 2 yi<i<n 

uniformly in the cylinder [0, T*) x Mo, where Mo is any manifold that can serve 
as a base manifold for the embedding x — x(t, £). 

The parabolic equations that are used for these curvature estimates are 
(|4.36| on pageHH usually for an upper estimate, and (|4.33p on page [14] for the 
lower estimate. Indeed, suppose that the flow starts at the upper barrier, then 

(6.27) F > f 
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at t = and this estimate remains valid throughout the evolution because of 
the parabolic maximum principle, use (|4.36j) . Then, if upper estimates for the 
Ki have been derived and if / > uniformly, then we conclude from (|6.27p that 
the Ki stay in a compact set inside the open cone r, since 

(6.28) F, or = 0. 

To obtain higher order estimates we are going to exploit the fact that the 
flow hypersurfaces are graphs over Sq in an essential way, namely, we look 
at the associated scalar flow equation (|4.21|) on page [12] satisfied by u. This 
equation is a nonlinear uniformly parabolic equation, where the operator <&(F) 
is also concave in hij, or equivalently, convex in Uij, i.e., the C 2, "-estimates of 
Krylov and Safonov, Chapter 5.5] or see [3S] Chapter 10.6] for a very clear 
and readable presentation, are applicable, yielding uniform estimates for the 
standard parabolic Holder semi-norm 

(6.29) [D 2 u]p.,Q T 
for some < j3 < a in the cylinder 

(6.30) Q=[0,T)xS , 

independent of < T < T* , which in turn will lead to H m+2+a - :s± ¥ L ^ (Q T ) 
estimates, cf. 18, Theorem 2.5.9, Remark 2.6.2]. 

H m + 2 + a ' m± 2 ±s - (Qt) is a parabolic Holder space, cf. [37J P- 7] for the original 
definition and [TS] Note 2.5.4] in the present context. 

The estimate (|6.29p combined with the uniform C 2 -norm leads to uniform 
C 2, ^(<5>o)-estimates independent of T. 

These estimates imply that T* = oo. 

Thus, it remains to prove that u(t, •) converges in C m+2 (iSo) to a stationary 
solution u, which is then also of class C m+2,a (So) in view of the Schauder 
theory. 

Because of the preceding a priori estimates u(t, •) is precompact in C 2 (Sq). 
Moreover, we deduce from the scalar flow equation (|4.21| on page [12] that u 
has a sign, i.e., the u(t,-) converge monotonely in C°(6>o) to u and therefore 
also in C 2 (S ). 

To prove that graph u is a solution, we again look at (|4.21[) and integrate it 
with respect to t to obtain for fixed x G So 

/oo 
e-*«|*-/|, 

where we used that (<P — f) has a sign, hence (<P — f)(t, x) has to vanish when 
t tends to infinity, at least for a subsequence, but this suffices to conclude that 
graph u is a stationary solution and 

(6.32) lim (<£-/) = 0. 

t— >oo 

Using the convergence of u to u in C 2 , we can then prove: 
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6.5. Theorem. The functions u(t, •) converge in C m+2 (<So) to u, if the data 
satisfy the assumptions in Theorem \6.2\ since we have 

(6.33) u^H m+2+ ^ s± ^ ±A {Q), 
where Q = Qao. 

Proof. Out of convention let us write a instead of [3 knowing that a is the 
Holder exponent in (|6.29p . 

We shall reduce the Schauder estimates to the standard Schauder estimates 
in for the heat equation with a right-hand side by using the already estab- 
lished results 1|6.29[) and 

(6.34) u(t,-) -> ueC m+2 ' a (S ). 

C 2 (So) 

Let (Uk) be a finite open covering of Sq such that each Uk is contained in a 
coordinate chart and 

(6.35) diamf/fc < p, 

p small, p will be specified in the proof, and let (i]k) be a subordinate finite 
partition of unity of class C m+2 ' a . 
Since 

(6.36) u e #™+2+«,^±i^ (Q T ) 
for any finite T, cf. jT5J Lemma 2.6.1], and hence 

(6.37) u(t, •) e C m+2 ' a {S a ) VO < t < oo 

we shall choose uq — u(to,-) as initial value for some large to such that 

(6.38) \a ij (t,-)-a^\ 0tSo < e/2 V* > t , 
where 

(6.39) a lj = v 2 &F ij 

and a lJ is defined correspondingly for M — graph u. 

However, making a variable transformation we shall always assume that 
to = and uq — u(0, •). 

We shall prove (|6.33p successively. 

(i) Let us first show that 

(6.40) D x ue H 2+a '^(Q). 

This will be achieved, if we show that for an arbitrary £ € C m+1 ' a (T 1 ' (5 )) 

(6.41) L P = D i ueH 2+a ^(Q) 1 

cf. EH Remark 2.5.11]. 

Differentiating the scalar flow equation (14.21(1 on page IT21 with respect to £ 
we obtain 

(6.42) ip-a ij cp ij + b i cp i + cip = f, 
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where of course the symbol / has a different meaning then in (14.211) . 

Later we want to apply the Schauder estimates for solutions of the heat flow 
equation with right-hand side. In order to use elementary potential estimates 
we have to cut off ip near the origin t = by considering 

(6.43) (p = tp9, 

where = 9{t) is smooth satisfying 



(6.44) 0(f) 



* > 1, 
t<\. 



This modification doesn't cause any problems, since we already have a priori 
estimates for finite t, and we are only concerned about the range 1 < t < oo. 
ip satisfies the same equation as ip only the right-hand side has the additional 
summand wO. 

Let rj = rjk be one of the members of the partition of unity and set 

(6.45) w — iprj, 

then w satisfies a similar equation with slightly different right-hand side 

(6.46) w - a lj w tj + b l Wi + cw = f 

but we shall have this in mind when applying the estimates. 

The w(t,-) have compact support in one of the LVs, hence we can replace 
the covariant derivatives of w by ordinary partial derivatives without changing 
the structure of the equation and the properties of the right-hand side, which 
still only depends linearly on ip and Dip. 

We want to apply the well-known estimates for the ordinary heat flow equa- 
tion 

(6.47) w- Aw = f 

where w is defined in ffi. x M. n . 

To reduce the problem to this special form, we pick an arbitrary xq € Uk, 
set zq = (Q,xq), z = (t,x) and consider instead of (|6.46p 

w — a 13 (zq)wh = f 

(6.48) V ^ 3 .. 

= [a lJ (z) - a 13 (z )]wij - b l Wi -cw + f, 

where we emphasize that the difference 

(6.49) \a^(z)-a»(z )\ 

can be made smaller than any given e > by choosing p = p(e) in (|6.35[) and 
to = to(e) in (|6.38f> accordingly. Notice also that this equation can be extended 
into R x R™, since all functions have support in {t > |}. 

Let < T < oo be arbitrary, then all terms belong to the required function 
spaces in Qt and there holds 

(6-50) [w] 2+a ,Q T < c[f] a ,Q T , 
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where c = c(n,a). The brackets indicate the standard unweighted parabolic 
semi-norms, cf. [18j Definition 2.5.2], which are identical to those defined in 
[2"T1 p. 7], but there the brackets are replaced by kets. 
Thus, we conclude 

[wh+a,Q T < c sup \a %] {z) - a v (z )\[D 2 w} at Q T + c[/] q ,q t 

(6.51) u k x(o,T) 

+ C 1 {[D 2 u] atQ+ + [Du] a ^Q T + [u] at Q T + Mo, Q T + I^Ho, Qt}> 

where c% is independent of T, but dependent on r/k- Here we also used the fact 
that the lower order coefficients and ip, Dp are uniformly bounded. 
Choosing now e > so small that 

(6.52) ce < \ 

and p, to accordingly such the difference in (|6.49j) is smaller than e, we deduce 
^ [w]2+a,Q T < 2c[/] Q: q t 

+ 2c 1 {[D 2 u] a . Q+ + [Du] a .Q T + [u] a .Q T + \w\ 0y Q T + \D 2 w\q. Qt }. 



Summing over the partition of unity and noting that £ is arbitrary we see 
that in the preceding inequality we can replace w by Du everywhere resulting 
in the estimate 

[Du] 2 + a ,Q T < Cl[f] a ,Q T 

(6.54) 

+ c 1 {[D 2 u} a ^Q T + [Du] a Q T + [u] a ^Q T + \Du\ , q t + \D 3 u\ 0: q t }, 

where c\ is a new constant still independent of T. 

Now the only critical terms on the right-hand side are |D 3 m|o,q t , which can 
be estimated by (|6.57p . and the Holder semi- norms with respect to t 

( 6 - 55 ) [Du]9L tt ,Q T + [u]f , t ,Q T . 

The second one is taken care of by the boundedness of ii, see (|4.2ip on page[T2l 
while the first one is estimated with the help of equation (|6.42p revealing 

(6.56) \Dii\ < c{sup|u| 3 , So + |/|o,Q T }, 

[0,T] 

since for fixed but arbitrary t we have 

(6.57) \u\3,s < e[£ ,3 w]a,5 +C e |lt|o,S , 

where c e is independent of t. 
Hence we conclude 

(6.58) \Du\ 2+a, Qt — Const 
uniformly in T. 

(ii) Repeating these estimates successively for 2 < I < m we obtain uniform 
estimates for 

m 

(6.59) 

1=2 
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which, when combined with the uniform C 2 -estimates, yields 

(6.60) \u(t, •)|m+2, Q ,5o < const 

uniformly in < t < oo. 

Looking at the equation (|4.21ll we then deduce 

(6.61) \u(t, Olm.a.So < const 
uniformly in t. 

(iii) To obtain the estimates for D\u up to the order 

(6.62) [ m + 2 2 +» ] 

we differentiate the scalar curvature equation with respect to t as often as 
necessary and also with respect to the mixed derivatives D r t D s x to estimate 

(6.63) ]T D\D%u 

l<2r+s<m+2+a 

using (|6.60[) , (|6.6ip and the results from the prior differentiations. 

Combined with the estimates for the heat equation in R x R™ these estimates 
will also yield the necessary a priori estimates for the Holder semi-norms in Q, 
where again the smallness of f|6.49|) has to be used repeatedly. □ 

6.6. Remark. The preceding regularity result is also valid in Riemannian 
manifolds, if the flow hypersurfaces can be written as graphs in a Gaussian 
coordinate system. In fact the proof is unaware of the nature of the ambient 
space. 

With the method described above the following existence results have been 
proved in globally hyperbolic spacetimes with a compact Cauchy hypersurface. 
Q C N is always a precompact domain the boundary of which is decomposed 
as in (|6.4p and (|6.5p into an upper and lower barrier for the pair (F, /). We 
also apply the stability results from Section [5] and the just proved regularity of 
the convergence and formulate the theorems accordingly. 

By convergence of the flow in C m+2 we mean convergence of the leaves 
M(t) = graph u(t, •) in this norm. 

6.7. Theorem. Let Mi, Mi be lower resp. upper barriers for the pair (H, f), 
where f G C m ' Q (J?) and the Mi are of class C m+2 ' a , 4 < m, < a < 1, then 
the curvature flow 

x=(H- f)v 

(6.64) V 

x(Q) = x Q , 

where xq is an embedding of the initial hypersurface Mo = M<x exists for all 
time and converges in C m+2 to a stable solution M of class C m+2,a of the 
equation 

(6.65) H {M = /, 

provided the initial hypersurface is not already a solution. 
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The existence result was proved in [TTJ Theorem 2.2], see also [TBI Theorem 
4.2.1] and the remarks following the theorem. Notice that / isn't supposed to 
satisfy any sign condition. 

For spacetimes that satisfy the timelike convergence condition and for func- 
tions / with special structural conditions existence results via a mean curvature 
flow were first proved in [§]. 

The Gaussian curvature or the curvature functions F belonging to the larger 
class (K*), see [TU] for a definition, require that the admissible hypersurfaces 
are strictly convex. 

Moreover, proving a priori estimates for the second fundamental form of 
a hypersurface M in general semi-Riemannian manifolds, when the curvature 
function is not the mean curvature, or does not behave similar to it, requires 
that a strictly convex function \ is defined in a neighbourhood of the hypersur- 
face, see Lemma 12.21 on page [3] where sufficient assumptions are stated which 
imply the existence of strictly convex functions. 

Furthermore, when we consider curvature functions of class (K*), notice 
that the Gaussian curvature belongs to that class, then the right-hand side / 
can be defined in T(f2) instead of fl, i.e., in a local trivialization of the tangent 
bundle / can be expressed as 

(6.66) f = f(x,u) A v G T X (N). 

We shall formulate the existence results with this more general assumption, 
though of course any stability claim only makes sense for / = f(x). 

6.8. Theorem. Let F £ C m ' a (r + ), 4 < m, < a < 1, be a curvature 
function of class (K*), let < / G C m ' Q (T(J?)), and let Mi, M 2 be lower resp. 
upper barriers for (F, f) of class C m+2 ' a . Then the curvature flow 

x = {$- f)v 
x(0) = x 

where 4>(r) = logr and xq is an embedding of Mq = Mi, exists for all time and 
converges in C rn+2 to a stationary solution M £ (jrn+2,a Q j ^ g e q Ua fa on 

(6.68) F lu = f 

provided the initial hypersurface M 2 is not already a stationary solution and 
there exists a strictly convex function \ G C 2 (Q). 

When f = f{x) and F is of class (D), then M is stable. 

The theorem was proved in 1 10] when / is only defined in Q and in the 
general case in [TH1 Theorem 4.1.1]. 

When F = H2 is the scalar curvature operator, then the requirement that 
/ is defined in the tangent bundle and not merely in N is a necessity, if the 
scalar curvature is to be prescribed. To prove existence results in this case, / 



(6.67) 
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has to satisfy some natural structural conditions, namely, 

(6.69) 0<ci<f(x,u) if (v,v)=-l, 

(6.70) lU{x,v)l <c 2 (l + |H| 2 ), 
and 

(6.71) |||/^(x^)||| <c 3 (l + |||HII), 

for all x G (2 and all past directed timelike vectors v G T x (f2), where ||| • ||| is a 
Riemannian reference metric. 

Applying a curvature flow to obtain stationary solutions requires to approx- 
imate F and / by functions F e and and to use these functions for the flow. 
The F e are the e-regularizations of F, which we already discussed before, cf. 
(15.301) on page [19] Let us also write F instead of F e as before. 

The functions have the property that |||/fc^||| only grows linearly in \\\v\\\ 
and IH/feyO (x, is bounded. To simplify the presentation we shall therefore 
assume that / satisfies 

(6.72) \\\f (x,v)\\\<c 2 (l + \\\v\\\), 

(6.73) \\\fuf{x,v)\\\ <c 3 , 
and also 

(6.74) 0<cx<f(x,v) G T X {N), (v,v) < 0, 

although the last assumption is only a minor point that can easily be dealt with, 
see [El Remark 2.6], and [HI Section 7 and 8] for the other approximations of 
/• 

The barriers Mi, i = 1,2, for (F, /) satisfy the barrier condition of course 
only weakly, i.e., no strict inequalities; however, because of the e-regularization 
we need strict inequalities, so that the Mi's are also barriers for (F, /), if e is 
small. In [13, Remark 2.4 and Lemma 2.5] it is shown that strict inequalities 
for the barriers may be assumed without loss of generality. 

Now, we can formulate the existence result for the scalar curvature operator 
F = H2 under these provisions. 

6.9. Theorem. Let f £ C" ha (T(f>)), 4 < m, < a < 1, satisfy the con- 
ditions (|6~72|) . ([673]) and (|6J4]) . and let M 1; M 2 be strict lower resp. upper 
barriers of class C m+2,a for (F, /). Let F be the e-regularization of the scalar 
curvature operator F, then the curvature flow for F 

x{0) = x 

where #(r) = and xo is an embedding of Mq = M 2> exists for all time and 
converges in C m+2 to a stationary solution M e G C m+2 ' a of 

(6-76) F, Me = / 



CURVATURE FLOWS IN SEMI-RIEMANNIAN MANIFOLDS 



37 



provided there exists a strictly convex function x € C 2 (f2) and < e is small. 
The M e then converge in C m+2 to a solution M £ C m+2 < a of 

(6-77) F lM = f. 

If f = f(x) and N Einstein, then M is stable. 

These statements, except for the stability and the convergence in C m+2 , are 
proved in [13] , 

6.10. Remark. Let us now discuss the pure mean curvature flow 
(6.78) x = Hv 

with initial spacelike hypersurface Mo of class C m+2,a , m > 4 and < a < 1. 
From the corresponding scalar curvature flow (|4.21|) on page[l2]we immediately 
infer that the flow moves into the past of Mo, if 

(6-79) iJ, Mo > 

and into its future, if 

(6.80) H luo < 0. 

Let us only consider the case ()6.79|) and also assume that Mq is not maximal. 
From the a priori estimates in [111 Section 3 and Section 4] we then deduce 
that the flow remains smooth as long as it stays in a compact set of N, and if 
a compact, spacelike hypersurface Mi of class C 2 satisfying 

(6.81) H lMi < 

lies in the past of Mo, then the flow will exist for all time and converge in 
(jm+2 |- a gta^g maximal hypersurface M, hence a neighbourhood of M can 
be foliated by CMC hypersurfaces, where those in the future of M have positive 
mean curvature, in view of Remark 15.71 on page 1241 

Thus, the flow will converge if and only if such a hypersurface Mi lies in the 
past of M). 

Conversely, if there exists a compact, spacelike hypersurface Mi in N sat- 
isfying (|6.81|) . and there is no stable maximal hypersurface in its future, then 
this is a strong indication that N has no future singularity, assuming that 
such a singularity would produce spacelike hypersurfaces with positive mean 
curvature. 

An example of such a spacetime is the (n + 1) -dimensional de Sitter space 
which is geodesically complete and has exactly one maximal hypersurface M 
which is also totally geodesic but not stable, and the future resp. past of M 
are foliated by coordinate slices with negative resp. positive mean curvature. 

To conclude this section let us show which spacelike hypersurfaces satisfy 
C 1 -estimates automatically. 
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6.11. Theorem. Let M = graphu| s be a compact, spacelike hypersurface 
represented in a Gaussian coordinate system with unilateral bounded principal 
curvatures, e.g., 

(6.82) Ki > kq Vi. 

Then, the quantity v = , 1 can be estimated by 

^/l — \Du\ 2 

(6.83) v < c(\u\,S ,aij,ip,K ), 

where we assumed that in the Gaussian coordinate system the ambient metric 
has the form as in (|6.ip . 

Proof. We suppose as usual that the Gaussian coordinate system is future 
oriented, and that the second fundamental form is evaluated with respect to 
the past directed normal. We observe that 



(6.84) \\Duf = g^muj 



-2*\D»\ 2 



= e-^ 



hence, it is equivalent to find an a priori estimate for ||-Du||. 
Let A be a real parameter to be specified later, and set 

(6.85) w = ilog||L»u|| 2 + Au. 

We may regard w as being defined on So; thus, there is xq £ Sq such that 

(6.86) w(xq) — supw, 



and we conclude 



1 



(6.87) = Wi = UijU 3 + Xut 

in xq, where the covariant derivatives are taken with respect to the induced 
metric <jry, and the indices are also raised with respect to that metric. 

Expressing the second fundamental form of a graph with the help of the 
Hessian of the function 

(6.88) e^v^hij = -u i3 - r °o w i u i ~ ^oi u j ~ ^oj u i ~ ifj- 
we deduce further 

A||Du|| 4 = -Uiju'u 3 

(6.89) = e-^vhijuV + r o ° ||L>u|| 4 

+ 2f °^|| J Du|| 2 +r? j u i u j . 

Now, there holds 

(6.90) u l = g ij Uj = e- 2 ^a ij Uj v- 2 , 
and by assumption, 

(6.91) hijtfv? > k q \\Du\\ 2 , 
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i.e., the critical terms on the right-hand side of (|6.89p are of fourth order in 
||Dm|| with bounded coefficients, and we conclude that ||-Dm|| can't be too large 
in Xq if we choose A such that 

(6.92) A < -c|||f - 1 

with a suitable constant c; w, or equivalently, ||-Du|| is therefore uniformly 
bounded from above. □ 

Especially for convex graphs over So the term v is uniformly bounded as 
long as they stay in a compact set. 

7. The inverse mean curvature flow 
Let us now consider the inverse mean curvature flow (IMCF) 

(7.1) x = -H- 1 p 

with initial hypersurface Mq in a globally hyperbolic spacetime N with compact 
Cauchy hypersurface Sq. 

N is supposed to satisfy the timelike convergence condition 

(7.2) R a pv a v P > V{v,v) = -1. 

Spacetimes with compact Cauchy hypersurface that satisfy the timelike con- 
vergence condition are also called cosmological spacetimes, a terminology due 
to Bartnik. 

In such spacetimes the inverse mean curvature flow will be smooth as long 
as it stays in a compact set, and, if H | > and if the flow exists for all time, 
it will necessarily run into the future singularity, since the mean curvature of 
the flow hypersurfaces will become unbounded and the flow will run into the 
future of Mo . Hence the claim follows from Remark 16.41 on page |2"51 

However, it might be that the flow will run into the singularity in finite 
time. To exclude this behaviour we introduced in [15] the so-called strong 
volume decay condition, cf . Definition 17.21 A strong volume decay condition is 
both necessary and sufficient in order that the IMCF exists for all time. 

7.1. Theorem. Let N be a cosmological spacetime with compact Cauchy hy- 
persurface Sq and with a future mean curvature barrier. Let Mo be a closed, 
connected, spacelike hypersurface with positive mean curvature and assume fur- 
thermore that N satisfies a future volume decay condition. Then the LMCF 
(17. 1|) with initial hypersurface Mo exists for all time and provides a foliation of 
the future D+(M Q ) of M . 

The evolution parameter t can be chosen as a new time function. The flow 
hypersurfaces M (t) are the slices {t = const} and their volume satisfies 

(7.3) \M(t)\ = \Mo\e-K 



Defining a new time function r by choosing 
(7.4) T = l-e""' 
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we obtain < r < 1, 

(7.5) |M(t)| = |Mo|(1-t)", 

and £/ie future singularity corresponds to r = 1 . 

Moreover, the length L(^) of any future directed curve 7 starting from M (r) 
is bounded from above by 

(7.6) L( 7 )<c(l-r), 

where c — c(n, Mo). Thus, the expression 1 — r caw 6e looked at as the radius 
of the slices {t = const} as w;eH as a measure of the remaining life span of the 
spacetime. 

Next we shall define the strong volume decay condition. 

7.2. Definition. Suppose there exists a time function a; such that the future 
end of N is determined by {tq < x° < b} and the coordinate slices M T — {x° — 
t} have positive mean curvature with respect to the past directed normal for 
tq < t < b. In addition the volume \M T \ should satisfy 

(7.7) lim|M T |=0. 

A decay like that is normally associated with a future singularity and we 
simply call it volume decay. If (gij) is the induced metric of M T and g — 
det(.gjj), then we have 

(7.8) log ff (r ,a;)-log ff (r,a;)= f 2e^H(s,x) VxeS , 

J T(, 

where H(r, x) is the mean curvature of M T in (t, x). This relation can be easily 
derived from the relation (|3.8j) on page[5]and Remark l3.6l on page[7l A detailed 
proof is given in |12j . 

In view of (|7.7p the left-hand side of this equation tends to infinity if t 
approaches b for a.e. x £ So, i.e., 

(7.9) lim / e" p H{s 1 x) = 00 for a.e. x € Sq. 

T - h JT0 

Assume now, there exists a continuous, positive function ip — ip(r) such that 

(7.10) e^r, x) > ?{t) V (r, x) G (r , b) x So, 
where 

(7.11) / ^(r)=oo, 



T 

then we say that the future of N satisfies a strong volume decay condition. 

7.3. Remark, (i) By approximation we may assume that the function ip 
above is smooth. 
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(ii) A similar definition holds for the past of N by simply reversing the time 
direction. Notice that in this case the mean curvature of the coordinate slices 
has to be negative. 

7.4. Lemma. Suppose that the future of N satisfies a strong volume decay 
condition, then there exist a time function x = x°(x°), where x° is the time 
function in the strong volume decay condition, such that the mean curvature H 
of the slices x = const satisfies the estimate 

(7.12) e*H>l. 

The factor e^ is now the conformal factor in the representation 

(7.13) ds 2 = e 2 ^(-(di ) 2 + a^dx'dx 1 ). 

The range of x is equal to the interval [0, oo), i.e., the singularity corre- 
sponds to x = oo. 

A proof is given in [151 Lemma 1.4]. 

7.5. Remark. Theorem 17. 1 1 c an be generalized to spacetimes satisfying 

(7.14) R a pv a v f3 > -A V>,z/>=-l 

with a constant A > 0, if the mean curvature of the initial hypersurface Mq is 
sufficiently large 




cf. [25j . In that thesis it is also shown that the future mean curvature barrier 
assumption can be dropped, i.e., the strong volume decay condition is sufficient 
to prove that the IMCF exists for all time and provides a foliation of the future 
of Mq. Hence, the strong volume decay condition already implies the existence 
of a future mean curvature barrier, since the leaves of the IMCF define such a 
barrier. 

8. The IMCF in ARW spaces 

In the present section we consider spacetimes N satisfying some structural 
conditions, which are still fairly general, and prove convergence results for the 
leaves of the IMCF. 

Moreover, we define a new spacetime ./V by switching the light cone and 
using reflection to define a new time function, such that the two spacetimes 
N and N can be pasted together to yield a smooth manifold having a metric 
singularity, which, when viewed from the region N is a big crunch, and when 
viewed from TV is a big bang. 

The inverse mean curvature flows in N resp. N correspond to each other via 
reflection. Furthermore, the properly rcscalcd flow in N has a natural smooth 
extension of class C 3 across the singularity into N. With respect to this natural 
diffeomorphism we speak of a transition from big crunch to big bang. 
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8.1. Definition. A globally hyperbolic spacetime N, dim N = n+ 1, is said 
to be asymptotically Roberts on- Walker (ARW) with respect to the future, if a 
future end of N, N + , can be written as a product N + — [a, b) x So, where iSo 
is a Riemannian space, and there exists a future directed time function t = x a 
such that the metric in N + can be written as 

(8.1) ds 2 = e 2 ^{-(dx ) 2 + a ij (x°,x)dx i dx j } ) 
where So corresponds to a; = a, ip is of the form 

(8.2) 4>(x°,x) = f(x°) + i/>(x° t x), 

and we assume that there exists a positive constant Co and a smooth Riemann- 
ian metric a%j on So such that 

(8.3) lim = Co A lim fy(T, x) = aij(x), 

r — yb t — yb 

and 

(8.4) lim /(t) = -oo. 

r— yb 

Without loss of generality we shall assume Co — 1. Then N is ARW with 
respect to the future, if the metric is close to the Robertson- Walker metric 

(8.5) ds 2 = e 2f {-dx° 2 + a ij (x)dx i dx j } 

near the singularity r = b. By close we mean that the derivatives of arbitrary 
order with respect to space and time of the conformal metric e~ 2 *g a p in (18. 1|) 
should converge to the corresponding derivatives of the conformal limit metric 
in (|8.5p when a: tends to b. We emphasize that in our terminology Robertson- 
Walker metric does not imply that (ct^-) is a metric of constant curvature, it is 
only the spatial metric of a warped product. 

We assume, furthermore, that / satisfies the following five conditions 

(8.6) -/'>0, 
there exists lu £ R such that 

(8.7) n + uj-2>0 A lim|/'| 2 e ( ™ + "- 2)/ = m > 0. 

r—>b 

Set 7 = |(n + lu — 2), then there exists the limit 

(8.8) hW/"+7|/'| 2 ) 

T^b 

and 

(8.9) IW+7l/f)l <c m \f'\ m Vm>l, 
as well as 

(8.10) \D?f\ < c m \f'\ m Vm>l. 

If So is compact, then we call N a normalized ARW spacetime, if 

(8.11) / v /det^~=|5 n |. 

J So 
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8.2. Remark, (i) If these assumptions are satisfied, then the range of t is 
finite, hence, we may — and shall — assume w.l.o.g. that 6 = 0, i.e., 

(8.12) a<r<0. 

(ii) Any ARW spacetime with compact iSo can be normalized as one easily 
checks. For normalized ARW spaces the constant m in (|8.7p is defined uniquely 
and can be identified with the mass of N, cf. 20J . 

(iii) In view of the assumptions on / the mean curvature of the coordinate 
slices M T — {x° = t} tends to oo, if t goes to zero. 

(iv) ARW spaces with compact So satisfy a strong volume decay condition, 
cf. Definition 17.21 on page [301 

(v) Similarly one can define N to be ARW with respect to the past. In this 
case the singularity would lie in the past, correspond to r = 0, and the mean 
curvature of the coordinate slices would tend to — oo. 



We assume that N satisfies the timelike convergence condition and that Sq 
is compact. Consider the future end N + of N and let Mq C N + be a spacelike 
hypersurface with positive mean curvature H\ Ma > with respect to the past 
directed normal vector v — it will become apparent in a moment why we use the 
symbols H and v and not the usual ones H and v. Then, as we have proved 
in the preceding section, the inverse mean curvature flow 

(8.13) x^-H^v 

with initial hypersurface Mq exists for all time, is smooth, and runs straight 
into the future singularity. 

If we express the flow hypersurfaces M(t) as graphs over So 

(8.14) M{t) = graph u(*,-)> 
then we have proved in [14] 

8.3. Theorem, (i) Let N satisfy the above assumptions, then the range of 
the time function x° is finite, i.e., we may assume that 6 = 0. Set 

(8.15) u = ue 7 *, 

where 7 = -^7, then there are positive constants C\,C2 such that 

(8.16) -c 2 <u<-ci<0, 

and u converges in C°°(So) to a smooth function, ift goes to infinity. We shall 
also denote the limit function by u. 

(ii) Let cjij be the induced metric of the leaves M(t), then the rescaled metric 

(8.17) e»% 
converges in C ca (So) to 

(8.18) (7m) y (-«) to-y. 
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(iii) The leaves M(t) get more umbilical, ift tends to infinity, namely, there 
holds 

(8.19) H^ 1 ^ — ^H5{\ < ce~ 27 *. 
In case n + u) — 4 > ; we even get a better estimate 

(8.20) \h{ - ±H6(\ < ce -^ (n+UJ - 4)t . 

To prove the convergence results for the inverse mean curvature flow, we con- 
sider the flow hypersurfaces to be embedded in N equipped with the conformal 
metric 

(8.21) ds 2 = -(dx ) 2 + a l] (x°,x)dx i dx j . 

Though, formally, we have a different ambient space we still denote it by the 
same symbol N and distinguish only the metrics g a p and g a p 

(8.22) g af} = e 2 ^g af3 

and the corresponding geometric quantities of the hypersurfaces , gtj , v resp. 
hij,gij,v, etc., i.e., the standard notations now apply to the case when N is 
equipped with the metric in (|8.2ip . 

The second fundamental forms hj and h\ are related by 

(8.23) e^hi=hi+$ a v a 6j 
and, if we define F by 

(8.24) F = e^H, 
then 

(8.25) F = H — nvf + nip a v a , 
where 

(8.26) v = v~\ 
and the evolution equation can be written as 

(8.27) x = -F~ x v, 
since 

(8.28) v = e-V 

The flow exists for all time and is smooth, due to the results in the preceding 
section. 

Next, we want to show how the metric, the second fundamental form, and 
the normal vector of the hypersurfaces M(t) evolve by adapting the general 
evolution equations in Section [3] on page |4] to the present situation. 
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(8.29) 




(8.30) 




and 




(8.31) 


K 


(8.32) 


hij 



8.4. Lemma. The metric, the normal vector, and the second fundamental 
form of M(t) satisfy the evolution equations 

Qij ^F h^, 
v = V M {-F- 1 )=gV{-F- 1 ) t x 3 , 

(-F- V )l + F^hfhl + F~ x R af)l8 v a x^xlg k i 
-- {-F-% - F^h^hkj + F- 1 R afhS v a xfvix s j . 

Since the initial hypersurface is a graph over iSo , we can write 

(8.33) M(t) = graphu(% o Vt G J, 

where u is defined in the cylinder R + x So. We then deduce from (18.271) . looking 
at the component a = 0, that u satisfies a parabolic equation of the form 

(8.34) u=j, 

where we emphasize that the time derivative is a total derivative, i.e. 

du 

(8.35) u= — +Uix\ 

at 

Since the past directed normal can be expressed as 

(8.36) (*/*) = -e-*v~ x {\,u% 
we conclude from (|8.34j) 

du v 

(8 - 37) dt = F- 

For this new curvature flow the necessary decay estimates and convergence 
results can be proved, which in turn can be immediately translated to corre- 
sponding convergence results for the original IMCF. 

Transition from big crunch to big bang 

With the help of the convergence results in Thcorem l8.31 we can rescale the 
IMCF such that it can be extended past the singularity in a natural way. 

We define a new spacetime N by reflection and time reversal such that the 
IMCF in the old spacetime transforms to an IMCF in the new one. 

By switching the light cone we obtain a new spacetime N. The flow equation 
in N is independent of the time orientation, and we can write it as 

(8.38) x = -U- X v = -(-H)~ x (-i>) = -H^v, 

where the normal vector v = — v is past directed in TV and the mean curvature 
H = —H negative. 
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Introducing a new time function x = ~x° and formally new coordinates 
(x a ) by setting 

(8.39) x° = -x°, x l =x\ 

we define a spacetime N having the same metric as N — only expressed in the 
new coordinate system — such that the flow equation has the form 

(8.40) x = -H'H, 
where M{t) — graph u(t), u — —u, and 

(8.41) {v a )=-ee-^(l,u i ) 
in the new coordinates, since 

(8.42) 
and 

(8.43) fr* = -v\ 

The singularity in x a = is now a past singularity, and can be referred to 
as a big bang singularity. 

The union NUN is a smooth manifold, topologically a product (— a, a) x So — 
we are well aware that formally the singularity {0} x Sq is not part of the union; 
equipped with the respective metrics and time orientation it is a spacetime 
which has a (metric) singularity in ir = 0. The time function 



144) x 



a 



in N, 
in N, 



is smooth across the singularity and future directed. 

NUN can be regarded as a cyclic universe with a contracting part N = 
{x° < 0} and an expanding part N = {x° > 0} which are joined at the 
singularity {x° — 0}. 

It turns out that the inverse mean curvature flow, properly rescaled, defines 
a natural C 3 - diffeomorphism across the singularity and with respect to this 
diffeomorphism we speak of a transition from big crunch to big bang. 

Using the time function in (|8.44[) the inverse mean curvature flows in N and 
iV can be uniformly expressed in the form 

(8.45) 'x = ~H- X v, 

where (|8.45|) represents the original flow in N, if x° < 0, and the flow in (|8.40|) . 
if x° > 0. 

Let us now introduce a new flow parameter 



.46) 



-7 1 e 7 *, for the flow in N, 
7 -1 e -7 *, for the flow in N, 
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and define the flow y — y(s) by y(s) = x(t). y — y(s,£) is then defined in 
[— 7~ 1 ,7~ 1 ] x So, smooth in {s ^ 0}, and satisfies the evolution equation 



■47) y' = £y 

In [14] we proved: 



s < 0, 
H^ua*, s > 0. 



8.5. Theorem. The flow y = y(s,£) is of class C 3 in (— 7 _1 ) x So and 
defines a natural diffeomorphism across the singularity. The flow parameter s 
can be used as a new time function. 

8.6. Remark. The regularity result for the transition flow is optimal, i.e., 
given any < a < 1, then there is an ARW space such that the transition flow 
is not of class C 3 ' a , cf. [19]. 

8.7. Remark. Since ARW spaces have a future mean curvature barrier, a 
future end can be foliated by CMC hypersurfaces the mean curvature of which 
can be used as a new time function., see [9] and [22]. In [8] we study this folia- 
tion a bit more closely and prove that, when writing the CMC hypersurfaces as 
graphs M T = graph ip(r, •) in the special coordinate system of the ARW space, 
where r is the mean curvature, of M T then 

(8.48) ri-tp) 1 ^ 1 -> const > 0, 
notice that <p < 0, and hence 

(8.49) lim^^T^ 1 Vx,yeS . 

r^oo ip( T , y) 

Moreover, the new time function 

(8.50) s = -t~ 9 , q 



1 + 7 

can be extended to the mirror universe N by odd reflection as a function of 
class C 3 across the singularity with non- vanishing gradient. 
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